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Abstract. We study heavy- ion fusion reactions at energies near the Coulomb barrier, in 
particular with neutron-rich radioactive ion beams. Dynamic microscopic calculations are 
carried out on a three-dimensional lattice using the Density-Constrained Time-Dependent 
Hartree-Fock (DC-TDHF) method. New results are presented for the 132 Sn+ 40 Ca system 
which are compared to 132 Sn+ 48 Ca studied earlier. Our theoretical fusion cross-sections agree 
surprisingly well with recent data measured at HRIBF. We also study the near- and sub-barrier 
fusion of 24,16 on 12 C which is important to determine the composition and heating of the 
crust of accreting neutron stars. 



1. Introduction 

The calculation of heavy-ion interaction potentials is of fundamental importance for the study of 
fusion reactions between stable and neutron-rich nuclei, and for the study of superheavy element 
production. We have developed a fully microscopic method to extract ion-ion potentials directly 
from the Time-Dependent Hartree-Fock (TDHF) time-evolution of the nuclear system. The only 
input is the effective NN interaction (Skyrme), and there are no adjustable parameters. 

Radioactive ion beam facilities have opened up the possibility to study fusion reactions of 
neutron-rich nuclei. In many cases, the dynamics of the neutron rich skin of these nuclei enhances 
the sub-barrier fusion cross section over that predicted by a simple static barrier penetration 
model, but in some cases suppression of fusion is also observed. Very recently, at the HRIBF 
facility a series of experiments has been carried out with radioactive 132 Sn beams on 0,48 Ca 
targets pp. Also, at the GANIL-SPIRAL facility a reaccelerated beam of 20 O was used to 
measure near- and sub-barrier fusion with a 12 C target [2]. 

The time-dependent Hartree-Fock (TDHF) theory provides a useful foundation for a fully 
microscopic many-body theory of large amplitude collective motion OH] including deep-inelastic 
and fusion reactions. Recently it has become feasible, for the first time, to perform TDHF 
calculations on a 3D Cartesian grid without any symmetry restrictions and with much more 
accurate numerical methods [21 El El [9] . At the same time the quality of effective interactions 
has also been substantially improved \10\ [TT| I12j . During the past several years, we have 
developed the DC-TDHF method for calculating heavy-ion potentials [13] . and we have applied 
this method to calculate fusion and capture cross sections above and below the barrier. So 
far, we have studied the systems 132 Sn+ 64 Ni [II], 64 Ni+ 64 Ni, 16 O+ 208 Pb fl5J, 132 ' 124 Sn+ 96 Zr, 
and we have studied the entrance channel dynamics of hot and cold fusion reactions leading 



to superheavy element Z = 112 [16J. Most recently, we have investigated sub-barrier fusion 
and pre-equilibrium giant resonance excitation between various tin + calcium isotopes [17] and 
calcium + calcium isotopes [18]. Last not least, we have studied sub-barrier fusion reactions 
between both stable and neutron-rich isotopes of oxygen and carbon [19] that occur in the 
neutron star crust. In all cases, we have found good agreement between the measured fusion 
cross sections and the DC-TDHF results. This is rather remarkable given the fact that the 
only input in DC-TDHF is the Skyrme effective N-N interaction, and there are no adjustable 
parameters. 

2. Unrestricted TDHF dynamics 

The TDHF equations of motion are obtained from the variational principle 



where H denotes the quantum many-body Hamiltonian of the system consisting of kinetic energy 
and Coulomb / nuclear two-body interactions. The main approximation in TDHF is that the 
many-body wave function <3?(i) is assumed to be a single time-dependent Slater determinant 
which consists of an anti-symmetrized product of single-particle wave functions 



The variational principle confined to the subspace of Slater determinants yields the TDHF 
equations for the single-particle wave functions 



where h = dE/dp denotes the mean- field Hamiltonian. 

In the present TDHF calculations we use the Skyrme SLy4 interaction [10] for the nucleons 
including all of the time-odd terms in the mean- field Hamiltonian [5] . The numerical calculations 
are carried out on a 3D Cartesian lattice. For 40 ' 48 Ca+ 132 Sn the lattice spans 50 fm along the 
collision axis and 30 — 42 fm in the other two directions, depending on the impact parameter. 
Derivative operators on the lattice are represented by the Basis-Spline collocation method. One 
of the major advantages of this method is that we may use a relatively large grid spacing of 
1.0 fm and nevertheless achieve high numerical accuracy. First we generate very accurate static 
HF wave functions for the two nuclei on the 3D grid. The static HF equations are solved with 
the damped gradient iteration method. The initial separation of the two nuclei is 22 fm for 
central collisions. In the second step, we apply a boost operator to the single-particle wave 
functions. The time-propagation is carried out using a Taylor series expansion (up to orders 
10 — 12) of the unitary mean-field propagator, with a time step At = 0.4 fm/c. 

In Figures 11121.3141 we show contour plots of the TDHF mass density for a central collision of 
48 Ca+ 132 Sn at E c , m , = 140 MeV. The density distributions are shown function of time, or 
equivalently, as a function of the internuclear distance R. 

3. DC-TDHF method: ion-ion potential and fusion cross section 

In the absence of a true quantum many-body theory of barrier tunneling, sub-barrier fusion 
calculations are reduced to the calculation of a potential barrier between the interacting nuclei 
and a subsequent calculation of tunneling through the barrier jT4]. To date, theoretical studies 
of fusion cross sections are still dominated by phenomenological methods such as the coupled- 
channels (CC) approach [20J. In this approach, one either uses empirical parameterizations of 
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Figure 1. Central collision of 48 Ca+ 132 Sn at 
E c .m. = 140 MeV. Shown is a contour plot 
of the TDHF mass density at internuclear 
distance R = 22.0 fm. 
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Figure 2. Mass density at R = 11.8 fm. A 

neck between the two fragments is starting to 
develop. 
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Figure 3. Mass density at R = 10.4 fm. 
The neck broadens resulting in increased mass 
transfer. 
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Figure 4. Mass density at R = 9.66 fm 
corresponding to the capture point. 



the ion-ion potential (e.g. Woods-Saxon) or one calculates the potential with the double-folding 
method using experimental or theoretical nuclear densities for projectile and target. In the latter 
case, one relies on the "frozen density" or "sudden" approximation in which the nuclear densities 
are unchanged during the computation of the ion-ion potential as a function of the internuclear 
distance. The frozen density approximation ignores dynamical effects such as neck formation 
during the nuclear overlap. It has been demonstrated that for deep sub-barrier energies the 
inner part of the potential barrier plays a very important role [21j . 

While phenomenological methods provide a useful starting point for the analysis of fusion 
data, it is desirable to use a microscopic many-body theory. We have developed a fully 
microscopic method to extract heavy-ion interaction potentials V{R) from the TDHF time- 
evolution of the dinuclear system. 

In our DC-TDHF approach, the time-evolution takes place with no restrictions. At certain 
times t or, equivalently, at certain internuclear distances R{t) the instantaneous TDHF density 

PtdhfM) =< $(i)|p|$(t) > (4) 

is used to perform a static Hartree-Fock energy minimization 

5 < $ p \H - f d 3 r A(r) p(r) |$ p >= (5) 



while constraining the proton and neutron densities to be equal to the instantaneous TDHF 
densities 

< $p|p|$ p >= PtdhfM) • (6) 

These equations determine the state vector This means we allow the single-particle wave 
functions to rearrange themselves in such a way that the total energy is minimized, subject to 
the TDHF density constraint. We have a self-organizing system which selects its path following 
the dynamics given by the microscopic TDHF equations. 

In a typical DC-TDHF run, we utilize a few thousand time steps, and the density constraint 
is applied every 20 time steps. We refer to the minimized energy as the "density constrained 
energy" E DC (R) 

£bc(fl)=<$ p |%> . (7) 

The ion-ion interaction potential V(R) is essentially the same as Edc(R), except that it is 
renormalized by subtracting the constant binding energies E\ 1 and E& 2 of the two individual 
nuclei 



V(R) = E DC (R)-E Al -E A2 



(8) 



The interaction potentials calculated with the DC-TDHF method incorporate all of the 
dynamical entrance channel effects such as neck formation, particle exchange, internal 
excitations, and deformation effects. While the outer part of the potential barrier is largely 
determined by the entrance channel properties, the inner part of the potential barrier is strongly 
sensitive to dynamical effects such as particle transfer and neck formation. 

Using TDHF dynamics, it is also possible to compute the corresponding coordinate dependent 
mass parameter M(R) using energy conservation at zero impact parameter [15]. As expected, 
at large distance R the mass M(R) is equal to the reduced mass \x of the system. At smaller 
distances, when the nuclei overlap, the mass parameter increases in all cases. Instead of solving 
the Schrbdinger equation with coordinate dependent mass parameter M{R) it is numerically 
advantageous to use the constant reduced mass \i and to transfer the coordinate-dependence of 
the mass to a scaled potential using a scale transformation (for details, see [15]). For simplicity 
of notation, we denote this transformed potential in the following by V(R). In general, we 
observe that the coordinate-dependent mass changes only the interior region of the potential 
barriers, and this change is most pronounced at low E c . m , energies. 

The Schrbdinger equation for the coordinate R involving constant reduced mass \i and 
transformed potential V(R) has the familiar form 



-h 2 d 2 h 2 £(£ + l) 
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+ V(R) - E c 
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(9) 



By numerical integration of Eq. ([9]) using the well-established Incoming Wave Boundary 
Condition (IWBC) method we obtain the barrier penetrabilities Tg which determine the total 
fusion cross section 

tt?i 2 00 

0-fuB^c.m.) = — E( 2 ^ + l)ZK^cm.) • (10) 

1=0 



4. Numerical results 

In Fig. [5]we show heavy-ion interaction potentials for 132 Sn+ 40 Ca calculated with the DC-TDHF 
method. The potentials have been transformed so that they correspond to a constant reduced 
mass [i. In general, the potentials are energy-dependent and they are shown here at four different 
TDHF energies. Our results demonstrate that in these heavy systems the potential barrier height 
increases dramatically with increasing energy -EtdhFi and the barrier peak moves inward towards 
smaller i?-values. Note that the potential calculated at high energy (-Etdhf = 180 MeV) has a 
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Figure 5. Transformed heavy- ion interaction 
potentials V(R) corresponding to the reduced 
mass \x. The potentials show a dramatic 
dependence on energy. 
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Figure 6. Total fusion cross sections 
obtained with the DC-TDHF method for 
132 Sn+ 40 Ca. The cross sections are calculated 
from Eq. fjXOj) using the transformed potential 
V{R) in Fig. [5] at four different energies. The 
experimental data are taken from Ref. pQ. 



barrier Eq = 115.3 MeV located at B, = 11.8 fm, whereas the potential calculated at low energy 
(-&TDHF = 115 MeV) has a barrier of only E B = 111.5 MeV located at R = 12.4 fm. 

In Fig. O we show total fusion cross sections obtained with the DC-TDHF method for 
132 Sn+ 40 Ca. Also shown are experimental data obtained at HRIBF [1J. The main point of this 
display is to demonstrate that the energy-dependence of the heavy-ion potential is crucial for 
an understanding of the strong fusion enhancement at subbarrier energies. At very high energy 
(-Etdhf = 180 MeV) the potential approaches the limit of the frozen density approximation: 
the collision is so fast that the nuclei have no time to rearrange their densities. We observe 
that the measured fusion cross sections at energies E c , m _ > 118 MeV are well-described by this 
high-energy potential. At low energies, however, the situation is quite different: Due to the slow 
motion of the nuclei, the nuclear densities have time to rearrange resulting in neck formation, 
surface vibrations, particle transfer, and the interior region of the heavy-ion potential is strongly 
modified. These effects (which are included in TDHF) lower the fusion barrier as described 
above, thus strongly enhancing the subbarrier fusion cross section. We observe that the data 
points measured at the two lowest energies E c _ m _ = 108.6 and 111 MeV can be understood by 
the heavy-ion potential calculated at -Etdhf = 115 MeV. 

In Fig. [8] we show total fusion cross sections for 132 Sn+ 48 Ca which contains 8 additional 
neutrons. In this case, we have interpolated the theoretical cross sections obtained with the 
energy-dependent DC-TDHF potentials [T7j. We can see that our theoretical cross sections agree 
remarkably well with the experimental data. If one compares the fusion cross sections for both 
systems at low energies, one finds the surprising result that fusion of 132 Sn with 40 Ca yields a 
larger cross section than with 48 Ca. For example, at E c m ^ = 110 MeV we find an experimental 
cross section of ~ 6 mb for 132 Sn+ 40 Ca as compared to 0.8 mb for the more neutron-rich 
system 132 Sn+ 48 Ca. This behavior can be understood by examining the corresponding heavy-ion 
potentials which are shown in Fig. [JJ Both potentials have been calculated at the same center- 
of-mass energy -Etdhf = 120 MeV. We observe that while the barrier heights and positions 
for both systems are approximately the same, the width of the DC-TDHF potential barrier 
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Figure 7. Transformed heavy-ion potentials 
V(R) at low energy ^tdhf = 120 MeV for 
the systems 132 Sn+ 40 ' 48 Ca. 
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Figure 8. Total fusion cross sections 
obtained with the DC-TDHF method for 
132 Sn+ 48 Ca. The experimental data are 
taken from Ref. pQ. 



for 3 Sn+ °Ca is substantially smaller than for 3 Sn+ Ca, resulting in enhanced sub-barrier 
fusion at low energy. 

Fusion of very neutron rich nuclei may be important to determine the composition and 
heating of the crust of accreting neutron stars [12]. In Fig. Owe show the DC-TDHF potential 
barriers for the C+O system. The higher barrier corresponds to the 12 C+ 16 system and has 
a peak energy of 7.77 MeV. The barrier for the 12 C+ 24 system occurs at a slightly larger R 
value with a barrier peak of 6.64 MeV. Figure [TOl shows the corresponding cross sections for the 
two reactions. Also shown are the experimental data from Ref. [23]. The DC-TDHF potential 
reproduces the experimental cross-sections quite well for the 12 C+ 16 system, and the cross 
section for the neutron rich 12 C+ 24 is predicted to be larger than that for 12 C + 16 0. 

We have compared our results to the phenomenological Sao Paulo barrier penetration model 
which calculates an effective potential by folding over static densities for the projectile and 
target. There are interesting differences between our calculations and the Sao Paulo model. For 
the asymmetric systems 12 C+ O or 16 0+ 24 we predict an order of magnitude larger cross 
section than for the Sao Paulo model. This is likely due to dynamical effects that change the 
nuclear densities during the collision process caused by the rearrangement of the single-particle 
wave functions. This enhancement of fusion cross sections of very neutron rich nuclei can be 
tested in the laboratory with radioactive beams. 

5. Summary and conclusions 

There is mounting evidence that the TDHF equations give a good dynamical description of 
the early stages of low-energy heavy-ion reactions. We have developed powerful methods for 
extracting more information from the TDHF dynamical evolution, e.g. the ion-ion interaction 
potential V(R), the coordinate-dependent mass parameter M(R), and the dynamical excitation 
energy E*(R). In general, V(R) depends on the center-of-mass energy and, in the case of 
deformed nuclei, on the Euler orientation angles. The only input is the Skyrme N-N interaction, 
and there are no adjustable parameters. We have calculated fusion cross sections below and 
above the barrier for 18 different heavy-ion systems so far. Overall, there is surprisingly good 
agreement with the experimental fusion data. In the future, we believe that an effort is needed 
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Figure 9. DC-TDHF heavy-ion potentials 
for the systems 12 C+ 16 - 24 0. 
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Figure 10. Total fusion cross sections versus 
center-of-mass energy for fusion of carbon 
with oxygen isotopes. The experimental data 
are from Ref. 123 . 



to incorporate deformation and scattering information into the Skyrme N-N parametrization. 
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